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We review recent advances in the theory of strong-interaction effects and final- 
state interactions in hadronic weak decays of heavy mesons. In the heavy-quark 
limit, the amplitudes for most nonleptonic, two-body B decays can be calculated 
from first principles and expressed in terms of semileptonic form factors and light- 
cone distribution amplitudes. We summarize the main features of this novel QCD 
factorization and discuss its phenomenological applications to B — » D^*'L decays 
(with L a light meson), and to the rare charmless decays B — » ttK and B — > tttt. 



1 Introduction 

The theoretical description of hadronic weak decays is difficult due to nonper- 
turbative strong-interaction dynamics. This affects the interpretation of data 
collected at the B factories and in many cases limits our ability to uncover 
the origin of CP violation and search for New Physics. The complexity of the 
problem is illustrated in the cartoon on the left-hand side of Fig. [l| 

It is well known how to control the effects of hard gluons with virtuality 
between the electroweak scale Mw and the scale tub characteristic to the de- 
cays of interest. They can be dealt with by constructing a low-energy effective 
weak Hamiltonian 

H ^ = % E A ? KM C *(Mw/v) OiQj) + h.c, (1) 

* i 

where \f KM are products of CKM matrix elements, Ci(M\y / /*) are calculable 
short-distance coefficients, and Oi(/i) are local operators renormalized at a 
scale /i = O(mfl). The challenge is to calculate the hadronic matrix elements 
of these operators with controlled theoretical uncertainties, using a systematic 
approximation scheme. 

Previous field-theoretic attempts to evaluate these matrix elements have 
employed dynamical schemes such as lattice field theory, QCD sum rules, or 
the hard-scattering approach. The first two have difficulties in accounting 
for final-state rescattering, which however is important for predicting direct 
CP asymmetries. The hard-scattering approach misses the leading soft con- 
tribution to the B — > meson transition form factors and thus falls short of 
reproducing the correct magnitude of the decay amplitudes. In view of these 
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Figure 1: Left: Strong-interaction effects in a hadronic weak decay. Right: QCD factoriza- 
tion in the heavy-quark limit. The second term is power suppressed for B — » Diz, but must 
be kept for decays with two light mesons in the final state, such as B — > ttK. Contributions 
not shown (such as weak annihilation graphs) are power suppressed. 

difficulties, most previous analyses of hadronic decays have employed phe- 
nomenological models such as "naive" or "generalized factorization" , in which 
the complicated matrix elements of four-quark operators in the effective weak 
Hamiltonian are replaced, in an ad hoc way, by products of current matrix 
elements. Corrections to this approximation are accounted for by introduc- 
ing a set of phenomenological parameters a,;. A different strategy is to clas- 
sify nonlcptonic decay amplitudes according to flavor topologies ("trees" and 
"penguins"), which can be decomposed into SU(3) or isospin amplitudes. This 
leads to relations between decay amplitudes in the flavor-symmetry limit. No 
attempt is made, however, to compute these amplitudes from first principles. 

2 QCD Factorization Formula 

Here we summarize recent progress in the theoretical understanding of non- 
lcptonic decay amplitudes in the heavy-quark limitu'EI The underlying idea is 
to exploit the presence of a large scale, i.e., the fact that nib ^> Aqcd- In 
order to disentangle the physics associated with these two scales, we factorize 
and compute hard contributions to the decay amplitudes arising from gluons 
with virtuality of order mb, and parameterize soft and collinear contributions. 
Considering the cartoon in Fig. |l|, we denote by Mi the meson that absorbs 
the spectator quark of the B meson, and by M2 the meson at the upper ver- 
tex, to which we refer as the "emission particle". We find that at leading 
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power in Aqcd / mi all long-distance contributions to the decay amplitudes 
can be factorized into semileptonic form factors and meson light-cone distri- 
bution amplitudes, which are much simpler quantities than the nonleptonic 
amplitudes themselves. A graphical representation of the resulting "factoriza- 
tion formula" is shown on the right-hand side in Fig. |l|. The physical picture 
underlying factorization is color transparency!™ If the emission particle is a 
light meson, its constituents carry large energy of order mj and are nearly 
collincar. Soft gluons coupling to this system see only its net zero color charge 
and hence decouple. Interactions with the color dipole of the small gg-pair are 
power suppressed in the heavy-quark limit. 

For B decays into final states containing a heavy charm meson and a light 
meson the factorization formula takes the form 

1 

{D^L-\ 0i \B d ) =J2Ff- DM f L JduTl(u)* L (u)+o(^} , (2) 
j o 

where Oi is an operator in the effective weak Hamiltonian ([!]), F^^ Di ' are 
transition form factors (evaluated at q 2 = m 2 L « 0), Jl and $l are the decay 
constant and leading-twist light-cone distribution amplitude of the light meson, 
and T~ are process-dependent hard-scattering kernels. For decays into final 
states containing two light mesons there is a second type of contribution to 
the factorization formula, which involves a hard interaction with the spectator 
quark in the B meson. This is shown by the second graph on the right-hand side 
m Fig. 0. Below we focus first on B — > D^L decays (with L a light meson), 
where this second term is power suppressed and caiL-be neglected. Decays 
into two light final-state mesons are more complicatedlM and will be discussed 
briefly in Sect. || A more detailed account of the conceptual foundations of the 
QCD factorization approach has been presented by Beneke at this Workshops. 

The factorization formula for nonleptonic decays provides a model-inde- 
pendent basis for the analysis of these processes in an expansion in powers 
and logarithms of Aqcd/?™;)- At leading power, but to all orders in a s , the 
decay amplitudes assume the factorized form shown in (|J). Having such a 
formalism based on power counting in Aqcd/^6 is of great importance to 
the theoretical description of hadronic weak decays, since it provides a well- 
defined limit of QCD in which these processes admit a rigorous, theoretical 
description. (For instance, the possibility to compute systematically 0(a s ) 
corrections to "naive factorization" , which emerges as the leading term in the 
heavy-quark limit, solves the old problem of renormalization-scale and scheme 
dependences of nonleptonic amplitudes.) The usefulness of this new scheme 
may be compared with the usefulness of the heavy-quark effective theory for 
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the analysis of exclusive semileptonic decays of heavy mesons, or of the heavy- 
quark expansion for the analysis of inclusive decay rates. In all three cases, it is 
the fact that hadronic uncertainties can be eliminated up to power corrections 
in Aqcd /nib that has advanced our ability to control theoretical errors. 

It must be stressed, however, that we are just beginning to explore the 
theory of nonleptonic B decays. Some important conceptual problems remain 
to be better understood. In the next few years it will be important to further 
develop the approach. This should include an all-orders proof of factoriza- 
tion at leading power, the development of a formalism for dealing with power 
corrections to factorization, understanding the light-cone structure of heavy 
mesons, and understanding the relevance (or irrelevance) of Sudakov form fac- 
tors. Also, we must gauge the accuracy of the approach by learning about the 
magnitude of corrections to the heavy-quark limit from extensive comparisons 
of theoretical predictions with data. 

As experience with previous heavy-quark expansions has shown, this is 
going to be a long route. Yet, already we have obtained important insights. 
Before turning to specific applications, let us mention three points here: 

1. Corrections to "naive factorization" (usually called "nonfactorizable 
effects") are process dependent, in contrast with a basic assumption underlying 
models of "generalized factorization" . 

2. The physics of nonleptonic decays is both rich and complicated. There 
may, in general, be an interplay of several small parameters (Wilson coeffi- 
cients, CKM factors, 1/N C , etc.) in addition to the small parameter Aqcd/^6 
relevant to QCD factorization. Also, several not-so-well-known input param- 
eters (e.g., heavy-to-light form factors and light-cone distribution amplitudes) 
introducing numerical uncertainties in the predictions. 

3. Strong-interaction phases arising from final-state interactions are sup- 
pressed in the heavy-quark limit. More precisely, the imaginary parts of non- 
leptonic decay amplitudes are suppressed by at least one power of a s (rrib) or 
Aqcd/^6- At leading power, the phases are calculable from the imaginary 
parts of the hard-scattering kernels in the factorization formula. 



3 Applications to B d —* D^+L Decays 

Our result for the nonleptonic Bd — * D^ + L~ decay amplitudes (with L a 
light meson) can be compactly expressed in terms of the matrix elements of a 
"transition operator" 



r = ^v: d v cb 



ai {DL)Q v ~ ai (D*L)Q A 



(3) 
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where the hadronic matrix elements of the operators Qv — c^b <g> dy^l — 75)1* 
and Qa — 27^756 <g> ^7^(1 — 75 )u are understood to be evaluated in factor- 
ized form. Eq. (||) defines the quantities a\{D^L), which include the leading 
"nonfactorizable" corrections, in a renormalization-group invariant way. To 
leading power in Aqcd /mb these quantities should not be interpreted as phe- 
nomenological parameters (as is usually done) , because they are dominated by 
hard gluon exchange and thus calculable in QCD. At next-to-leading order in 
a s we obtain □ 



«i 



(D«L) = Ci(m b ) + 



C 2 (m b ) 



C F a s (m b ) 

47T 



1 

J duF(u,±z)<P L (u) 



, (4) 



where Ci(rrn,) are the so-called "renormalization-scheme independent" Wilson 
coefficients, and the upper (lower) sign in the second argument of the function 
F(u,±z) refers to a D (D*) meson in the final state. The exact expression 
for this function is known but not relevant to our discussion here. Note that 
the coefficients ai(DL) and a\{D*L) are nonuniversal, i.e., they are explicitly 
dependent on the nature of the final-state mesons. Politzer and Wise have 
computed the "nonfactorizable" vertex corrections to the decay rate ratio of 
B d -> D+tt~ and B d -> D* + 7r-|3 This requires the symmetric part (with 
respect to u <-> 1 — u) of the difference F(u, z) — F(u, —z). We agree with their 
result. 

The expressions for the decay amplitudes obtained by evaluating the had- 
ronic matrix elements of the transition operator T involve products of CKM 
matrix elements, light-meson decay constants, B — > D^*' transition form fac- 
tors, and the QCD parameters a\{D^*'L). A numerical analysis shows that 
1 0,1 j = 1.055 ±0.025 for the decays considered below. Below we will use this as 
our central value. 



3.1 Tests of factorization 

A particularly clean test of our predictions is obtained by relating the B d — * 
D* + L~ decay rates to the differential semileptonic Bd — > D* + l~v decay rate 
evaluated at g? — m 2 L . In this way the parameters |ai| can be measured 
directly, since □ 

T(B d -» D*+L~) _ a _ 2nr a a 



dT(B d -> D*+l-D)/dq 2 



^\V ud \*fi\ ai {D*L)\*. (5) 



With our result for ai this relation becomes a prediction based on first prin- 
ciples of QCD. This is to be contrasted with the usual interpretation of this 
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formula, where a\ plays the role of a phenomenological parameter that is fitted 
from data. n 

Using results reported by the CLEO Collaboration^ we find 

K(£>*7r)| = 1.08±0.07 , \ai{D*p)\ = 1.09±0.10 , \ay{D*a{)\ = 1.08±0.11 , 

(6) 

in good agreement with our prediction. It is reassuring that the data show 
no evidence for large power corrections to our results. However, a further im- 
provement in the experimental accuracy would be desirable in order to become 
sensitive to process-dependent, nonfactorizable effects. 

3.2 Predictions for class-I decay amplitudes 

We now consider a larger set of so-called class-I decays of the form Bd — ► 
DW+Z,-, all of which are governed by the transition operator ([?]). In Tab. [I] 
we compare the QCD factorization predictions with experimental data. As 
previously we work in the heavy-quark limit, i.e., our predictions are model 
independent up to corrections suppressed by at least one power of Aqcd/?™&- 
There is good agreement between our predictions and the data within exper- 
imental errors, which however are still large. It would be desirable to reduce 
these errors to the percent level. (Note that we have not attempted to adjust 
the semileptonic form factors F+^ D and Aq^ d entering our results so as to 
obtain a best fit to the data.) 

The observation that the experimental data on class-I decays into heavy- 
light final states show good agreement with our predictions may be taken as 
(circumstantial) evidence that in these, decays there are no unexpectedly large 
power corrections. In our recent workEI we have addressed the important ques- 
tion of power corrections theoretically by providing estimates for two sources 
of power-suppressed effects: weak annihilation and spectator interactions. We 
stress that a complete account of power corrections to the heavy-quark limit 
cannot be performed in a systematic way, since these effects are no longer 
dominated by hard gluon exchange. However, we believe that our estimates 
are nevertheless instructive. 

We parameterize the annihilation contribution to the Bd — ► D + ir~ decay 
amplitude in terms of an amplitude A such that A{Bd — > D + tt~) = T + 
A, where T is the "tree topology" , which contains the dominant factorizable 
contribution. We find that A/T ~ 0.04. We have also obtained an estimate 
of nonfactorizable spectator interactions, which are part of T, finding that 
?spec/Ziead ~ —0.03. In both cases, the results exhibit the expected linear 
power suppression in the heavy-quark limit. We conclude that the typical size 
of power corrections in class-I decays into heavy-light final states is at the level 
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Table 1: Model-independent predictions for the branching ratios (in units of 10 — 3 ) of 
— > Z)(*' + L~ decays in the heavy-quark limit. Predictions are in units of ( | V c t, | /0.04) 2 X 
(|ai |/1.()5) 2 X (rg /1.56 ps). We showr-ejxperimental results reported by the CLEO Collab- 
oration □ and the Particle Data Grouptj 



Decay mode 


Theory (HQL) 


CLEO data PDG98 


B d -> D+7T- 
B d -» D+K- 
B d D+p- 
B d -» D + K*~ 
B d .D+a- 


3.27 x [F+(0)/0.6] 2 
0.25 x [F+(0)/0.6] 2 
7.64 x [F+(0)/0.6] 2 
0.39 x [F+(0)/0.6] 2 
7.76 x [F+(0)/0.6] 2 


2.50 ±0.40 3.0 ±0.4 
7.89 ± 1.39 7.9 ±1.4 
8.34 ± 1.66 6.0 ±3.3 


B d -» #*+tt- 
B d -> D*+K~ 
B d 

S d -> D*+K*- 
B d -> J D*+a- 


3.05 x L4 (0)/0.6] 2 
0.22 x [A o (0)/0.6] 2 
7.59 x [A o (0)/0.6] 2 
0.40 x [A o (0)/0.6] 2 
8.53 x [A o (0)/0.6] 2 


2.34 ± 0.32 2.8 ±0.2 
7.34 ± 1.00 6.7 ±3.3 
11.57 ±2.02 13.0 ±2.7 



of 10% or less, and thus our predictions for the values and the near universality 
of the parameters ai governing these decay modes appear robust. 

4 QCD Factorization in Charmless Hadronic B Decays 

The observation of rare charmless B decays into ttK and mr final states has 
resulted in a large amount of theoretical and phenomenological work that at- 
tempts to interpret these observations in terms of the factorization approxi- 
mation, or in terms of general parameterizations of the decay amplitudes. A 
detailed understanding of these amplitudes would help us to pin down the 
value of the CKM angle 7 = &rg(V* b ) using only data on CP-averaged branch- 
ing fractions. Here we briefly summarize the most important consequences of 
the QCD factorization approach for the nK and nn final statesB 

To leading order in an expansion in powers of Aqc D / 'm& , the B — » nK 
matrix elements obey the factorization formula shown on the right-hand side 
in Fig. 0: 

(ttK\ O t \B) = F^f K T l K>i * <S> K + F*-+ K U Tl,i * <&. 

+ fsfxU Tj l * $ B * $ K * , (7) 

where the *-products imply an integration over the light-cone momentum frac- 
tions of the constituent quarks inside the mesons. Our results are based on 
hard-scattering kernels including all corrections of order a s . Compared to our 
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previous discussion of B — ► irir decaysP the present analysis incorporates three 
new ingredients: the matrix elements of electroweak penguin operators (for 
ttK modes), hard-scattering kernels for asymmetric light-cone distributions, 
and the complete set of "chirally enhanced" 1 /rrib corrections. rTjhfi second ^nd 
third items have not been considered in other generalizations EJO of Ref. El to 
the ttK final states. The third one, in particular, is essential for estimating 
some of the theoretical uncertainties of the approach. For completeness, we 
note that the predictions from QCD factorization differ in essential aspects 
from those obtainecLin the conventional hard-scattering approachEj 

Following Ref. El, we have obtained the coefficients ai(nK) (with i — 
1, . . . , 10) of the effective, factorized "transition operator" defined analogously 
to the case of B — ► ittt decays, but augmented by coefficients related to elec- 
troweak penguin contributions. Chirally enhanced corrections arise from twist- 
3 two-particle light-cone distribution amplitudes, whose normalization involves 
the quark condensate. The relevant parameter, 2/j^/nib = —4(qq)/(f„mb), is 
formally of order Aqcd/w;,, but large numerically. The coefficients and as 
are multiplied by this parameter. There are also additional chirally enhanced 
corrections to the spectator-interaction term in (Q), which turn out to be the 
more important effect. In both cases, these corrections involve logarithmically 
divergent integrals, which violate factorization. For instance, for matrix el- 
ements of V — A operators the hard spectator interaction is proportional to 
(u = 1 — u, v = 1 — v) 

l 

dudv ^ , . / . . 2uir u \ , . 

— —$*-(«) **(») + —- ( 8 ) 
u v \ nib u J 

o 

when the spectator quark goes to the pion. (Here we used that the twist- 
3 distribution amplitudes can be taken to be the asymptotic ones when one 
neglects twist-3 corrections without the chiral enhancement.) The divergence 
of the u-integral in the second term as v — > implies that it is dominated by 
soft gluon exchange between the spectator quark and the quarks that form the 
kaon. We therefore treat the divergent integral X — f (dv/v) as an unknown 
parameter (different for the penguin and hard-scattering contributions), which 
may in principle be complex owing to soft rescattering in higher orders. In 
our numerical analysis we set X = ln(ms/0.35 GeV) + r, where r is chosen 
randomly inside a circle in the complex plane of radius 3 ("realistic" Lor 6 
("conservative"). Our results also depend on the i?-meson parameter El A^, 
which we vary between 0.2 and 0.5 GeV. Finally, there is in some cases a 
nonnegligible dependence of the coefficients ai(irK) on the renormalization 
scale, which we vary between mt/2 and 2m&. 



We take \V ub /V cb \ = 0.085 and m s (2GeV) = 110 MeV as fixed input to 
our analysis, noting that ultimately these Standard Model parameters, along 
with the CP-violating phase 7, might be extracted from a simultaneous fit to 
the B — > irK and B — > irir decay rates. We now summarize our main results. 



4-1 Results on SC/(3) breaking 

Bounds on 7 derived from ratios of CP-averaged B — » irK decay rates^Sli^ 
as well as the determination of 7 using the method of Ref. 113, rely on an 
estimate of ST/ (3) flavor-symmetry violations. We find that "nonfactorizable" 
St/(3)-breaking effects (i.e., effects not accounted for by the different decay 
constants and form factors of pions and kaons in the conventional factorization 
approximation) do not exceed the level of a few percent. 



4-2 Amplitude parameters 

The approach discussed here allows us to calculate the amplitudes for B — > 
ttK and B — > ttt: decays in terms of form factors and light-cone distribution 
amplitudes. Here we focus on decays whose branching ratios have already been 
measured. We write 



A(B° TT+Tr") = T [ e * + (P/TU] 
and parameterize the B — > irK amplitudes by0 



(9) 



A(B+ - 


-+ TT+K°) 


= P 


(1 - e a e^e*") 






V2A(B+ - 


+ TT Q K+) 


= P 




„ i<b f Art ii>j\ 

£3/2 e' F {e ' -qe ) 




-A{B° - 


> TT~K + ) 


= P 


1 - e a e l "e 17 - 




V2A{B° - 


-+ ir°K°) 


= A{B+ -> ir+K°) 












f ). 





(10) 

Table ^ summarizes the numerical values for the amplitude parameters for the 
conservative variation of X , and variation of the other parameters as explained 
above. The leading-order results correspond to the conventional factorization 
approximation at the fixed scale fi — m b . They are strongly scale dependent. 
In comparison, the scale-dependence of the next-to-leading order results is 
small, with the exception of qc e*" c . We stress that the ranges shown may 
overestimate the theoretical uncertainty, since the parameter X may ultimately 
be constrained from a subset of branching fractions. This is true, hp-particular, 
for the quantity e 3 / 2 in Tab. ||, which can be extracted from datat3 
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plble 2: Parameters for the B — ♦ 7T7r and B — » ^Js" decay amplitudes as denned in (|Sj) and 
hot), for conservative variation of all input parameters (see text). 



Parameter 


Range, NLO 


LO 


-e a e ir > 


(0.017-0.020) e l [ 13 ' 21 l° 


0.02 


£3/2 e** 


(0.20-0.38) e *H"V] 


0.36 




(0.53-0.63) e 1 ^ 7 ' 3 ! 


0.64 




(0.20-0.29) e 4 !- 19 - 3 ! 


0.33 




(0.00-0.22) e * [-180,180]° 


0.06 


(-P/2V 


(0.19-0.29) e 4 1- 1 ' 23 ! 


0.16 



4-3 Ratios of CP-averaged rates 

Since the relevant form factor F^~*^(0) is not well known, we consider only 
ratios of CP-averaged branching ratios here. We display these as functions 
of the weak phase 7 in Fig. ^|. For comparison, we show the data on the 
various ratios obtained using results reported by the CLEO CollaborationEl 
We also indicate the very recent results reported by the BaBar and Belle Col- 
laborations at the ICHEP2000 ConferenceLL3 From the plots and Tab. | it is 
apparent that the corrections with respect to the conventional factorization 
approximation are significant (and important to reduce the renormalization- 
scale dependence). Despite this fact, the qualitative pattern that emerges for 
the set of B — ► irK and 7T7r decay modes is similar to that in conventional 
factorization. In particular, the penguin-tree interference is constructive (de- 
structive) in B — > 7r + 7r~ (B — > ir~K + ) decays if 7 < 90°. Taking the currently 
favored range 7 = (60 ± 20)°, we find the following robust predictions: 



Br(7r+7r-) 
Br( 7 rT J FC ± ) 

B^TT^i) 

2Br(7r°if°) 
2Bt{tt°K ± ) 
Br{ir±K°) 

T B + B^TT^i) 



0.5-1.9 [0.25 ±0.10], 

0.9-1.4 [0.59 ±0.27], 

0.9-1.3 [1.27 ±0.47], 

0.6-1.0 [1.00 ±0.30]. (11) 
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25 50 75 100 125 150 175 25 50 75 100 125 150 175 

y/deg y/deg 



Figure 2: Ratios of CP-averaged B — » 7r_ft" and 7r7r decay rates. The scattered points cover 
a realistic (dark) and conservative (light) variation of input parameters. The dashed curve 
is the result obtained using "naive factorization" . 



The first ratio is clearly in disagreement with current CLEO dataO shown in. 
square brackets. However, there is good agreement with the recent results £3 
reported by BaBar (0.74 ± 0.29) and Belle (0.36 ± 0.26). 

The near equality of the second and third ratios in ( pd| ) is a consequence of 
isospin symmetrvB We find Br(B -> ir°K°) = (4.5±2.5)xl0- 6 x(|Kb|/0-04) 2 x 
[i 7 '^ 3 ^ 7r (0)/0.3] 2 almost independents of 7. This is three time smaller than 
the central value reported by CLEOEZI, (14.6 ±1;?±1^ x 10" 6 , and four times 
smaller than the central value reported by BELLE M, (21.0 ^Jt^H) x 10 ~ 6 - 
It will be interesting to follow how the comparison between data and theory 
will develop as the data become more precise. 

4-4 CP asymmetry in B — » 7r + 7r~ decay 

The stability of the prediction for the B — » tt + ii~ amplitude suggests that 
the CKM angle a can be extracted from the time-dependent mixing-induced 
CP asymmetry in this decay mode, without using isospin analysis. Figure |^ 
displays the coefficient S of — sin(Arns <J t) as a function of sin(2a) for sin(2/3) = 
0.75. For some values of S there is a two-fold ambiguity (assuming all angles are 
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sin(2(3)=0.75 



1 

0.75 
0.5 
0.25 


-0.25 
-0.5 
-0.75 





OC=45 



oc=135 



-0.75-0.5-0.25 0.25 0.5 0.75 

sin(2a0 



Figure 3: Mixing-induced CP asymmetry in B — > 7r+7r decays. The lower band refers to 
values 45° < a < 135°, the upper one to a < 45° (right) or a > 135° (left). We assume 
a, ft 7 e [0, 180°]. 

between 0° and 180°). A consistency check of the approach could be obtained, 
in principle, from the coefficient of the cos(Ams d t) term, which is given by 
the direct CP asymmetry in this decay. 

5 Summary and Outlook 

With the recent commissioning of the B factories and the planned emphasis 
on heavy- flavor physics in future collider experiments, the role of B decays in 
providing fundamental tests of the Standard Model and potential signatures 
of New Physics will continue to grow. In many cases the principal source of 
systematic uncertainty is a theoretical one, namely our inability to quantify 
the nonperturbative QCD effects present in these decays. This is true, in 
particular, for almost all measurements of direct CP violation. Our work 
provides a rigorous framework for the evaluation of strong-interaction effects 
for a large class of exclusive, two-body nonleptonic decays of B mesons. It gives 
a well-founded field-theoretic basis for phenomenological studies of exclusive 
hadronic B decays and a formal justification for the ideas of factorization. 

It is our belief that the factorization formula (j^) and its generalization to 
decays into two light mesons will form a useful basis for future phenomenolog- 
ical studies of nonleptonic B decays. We stress, however, that a considerable 
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amount of conceptual work remains to be completed. Theoretical investiga- 
tions along the lines discussed here should be pursued with vigor. We are 
confident that, ultimately, this research will result in a theory of nonleptonic 
B decays. 
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